Energy Gaps in Graphene Nanoribbons 
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Based on a first-principles approach, we present scaUng rules for the band gaps of graphene 
nanoribbons (GNRs) as a function of their widths. The GNRs considered have either armchair or 
zigzag shaped edges on both sides with hydrogen passivation. Both varieties of ribbons are shown 
to have band gaps. This differs from the results of simple tight-binding calculations or solutions of 
the Dirac's equation based on them. Our ab initio calculations show that the origin of energy gaps 
for GNRs with armchair shaped edges arises from both quantum confinement and the crucial effect 
of the edges. For GNRs with zigzag shaped edges, gaps appear because of a staggered sublattice 
potential on the hexagonal lattice due to edge magnetization. The rich gap structure for ribbons 
with armchair shaped edges is further obtained analytically including edge effects. These results 
reproduce our ab initio calculation results very well. 

PACS numbers: 73.22.-f, 75.70.Ak, 72.80.Rj 



The electronic structure of nanoscale carbon materi- 
als such as fullerenes and carbon nanotubes has been 
the subject of intensive research during the past two 
decades [ij because of fundamental scientific interest 
in nanomaterials and because of their versatile elec- 
tronic properties that are expected to be important 
for future nanoelectronics [1, Among the carbon 
nanostructures, a simple variation of graphene, ribbons 
with nanometer sized widths, has been studied exten- 
sively 0, H, @, S i, i, El El, [H, [3, [H. Because 
of recent progress in preparing single graphite layers on 
conventional device setups, the graphene nanoribbons 
(GNRs) with varying widths can be realized either by 
cutting mechanically exfoliated graphenes [l^ or 
by patterning epitaxially grown graphenes 1^, 2o| . 

Since GNRs are just geometrically terminated sin- 
gle graphite layers, their electronic structures have been 
modeled by imposing appropriate boundary conditions 
on Schrodinger's equation with simple tight-binding (TB) 
approximations based on 7r-states of carbon 0, H S 13 
or on a 2-dimensional free massless particle Dirac's equa- 
tion with an effective speed of Ught ('^ lO^m/s) [sHolllot. 
Within these models, it is predicted that GNRs with arm- 
chair shaped edges can be either metallic or semiconduct- 
ing depending on their widths ^, @, 0, S, H [lO] , and 
that GNRs with zigzag shaped edges are metallic with 
peculiar edge states on both sides of the ribbon regard- 
less of its widths 0, i, S, 0, i, S E E [II [11 ■ 



Although the aformentioned models are known to 
describe the low energy properties of graphene very 



well [H [11 [H, [13, [H 



|22|, a careful consideration of 
edge effects in nanometer sized ribbons are required to 
determine their bandgaps because, unlike the situation 
in graphene, the bonding characteristics between atoms 
change abruptly at the edges [3, [l3|. Moreover, the spin 
degree of freedom is also important because the GNRs 
with zigzag shaped edges have narrow-band edge states 



at the Fermi ener 
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implying possible magnetiza- 
Motivated by the recent 



experimental progress in this area, we have carried out 
first-principles calculation and theoretical analysis to ex- 
plore the relation between the bandgap and the geome- 
tries of GNRs. 

In this Letter, we show that GNRs with hydrogen 
passivated armchair or zigzag shaped edges both always 
have nonzero and direct bandgaps. The origins of the 
bandgaps for the different types of homogenous edges 
vary. The bandgaps of GNRs with armchair shaped edges 
originate from quantum confinement, and edge effects 
play a crucial role. For GNRs with zigzag shaped edges, 
the bandgaps arise from a staggered sublattice potential 
due to spin ordered states at the edges [1, 11, [13, [l3|- 
Although the ribbon widths and energy bandgaps of the 
GNRs are related to each other primarily in inverse pro- 
portion, there is a rich structure in the ratio of the pro- 
portionalities as in the behavior of carbon nanotubes [l[ . 
For GNRs with armchair edges, analytic scaling rules for 
the size of the bandgaps are obtained as a function of 
width including the effect of the edges and give a good 
agreement with our first-principles calculations. 

Our electronic structure calculation employs the first- 
principles self-consistent pseudopotential method [5^ us- 
ing the local (spin) density approximation (L(S)DA) [S^, 
241 . An energy cutoff (for a real space mesh size) of 400 
Rydbergs is employed and a double-*^ plus polarization 
basis set [2^ is used for the localized basis orbitals to 
deal with the many atoms in a unit cell of the GNRs 
of various widths. We obtained the electron density by 
integrating the density matrix with a Fermi-Dirac distri- 
bution (23^ [2^ . The geometry for each GNR studied is 
fully relaxed until the force on each of the atoms is less 
than 16pN. A k-point sampling of 32 (96) k-points that 
are uniformly positioned along the ID Brillouin zone is 
employed for GNRs with armchair (zigzag) shaped edges. 



FIG. 1: (a) Schematic of a 11-AGNR. The empty circles 
denote hydrogen atoms passivating the edge carbon atoms, 
and the black and grey rectangles represent atomic sites be- 
longing to different sublattice in the graphene structure. The 
ID unit cell distance and ribbon width are represented by da 
and Wa respectively. The carbon-carbon distances on the n- 
th dimer line is denoted by a„. (b) Schematic of a 6-ZGNR. 
The empty circles and rectangles follow the same convention 
described in (a). The ID unit cell distance and the ribbon 
width are denoted by dz and Wz respectively. 

We have tested the change of gap size by increasing the 
vacuum between edges from 20 to 40 A and between 
plances from 16 to 25 A and found no changes. 
Following previous convention 
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12l . Il3l . Il4l . |15|, the GNRs with armchair shaped edges 
on both sides are classified by the number of dimer lines 
(A^q) across the ribbon width as shown in Fig. 1(a). Like- 
wise, ribbons with zigzag shaped edges on both sides are 
classified by the number of the zigzag chains [Nz) across 
the ribbon width [Fig. 1(b)]. We refer to a CNR with Na 
dimer lines as a iVa-AGNR and a GNR with zigzag 
chains as a A^^-ZGNR. 

Our calculations show that the A^^-AGNRs are semi- 
conductors with energy gaps which decrease as a func- 
tion of increasing ribbon widths {wa). The variations 
in energy gap however exhibit three distinct family be- 
haviors [Fig. 2]. Moreover, the energy gaps obtained 
by a simple TB model are quite different from those 
by first-principles calculations. The TB results using 
a constant nearest neighbor hopping integral, t = 2.7 
eV 22| between 7r-electrons are summarized as function 
of width in Fig. 2(a). It shows that a A^a-AGNR is 
metallic if Na — ip + 2 (where p is a positive integer) or 
otherwise, it is semiconducting, in agreement with pre- 
vious calculations 0, d, H, 0,11, H, The gap of a 
A^a-AGNR (Atv^) is inversely proportional to its width, 
separated into basically two groups with a hierarchy of 
gap size given by > Asp+i > A3p+2(= 0) for aU 
p's. For the first-principles calculations, however, there 
are no metallic nanoribbons. The gaps as a function of 
ribbon width are now well separated into three differ- 
ent categories (or family structures) as shown in Fig. 
2(b). Moreover, the gap size hierarchy is also changed 
to Asp+i > Asp > A3p+2(7^ 0). For example, in the 
first-principles calculation for p = 13, the lowest energy 
gap is A38 — 45 meV and A40 — A39 = 68 meV, aU 
of which are quite larger values compared to those (0 
and —2 meV respectively) obtained from TB approxima- 
tions. The first-principles band structures of JV^-AGNRs 
are shown in Fig. 2(c) for A'a=12, 13, and 14. They 




FIG. 2: (color online) The variation of bandgaps of Na- 
AGNRs as a function of width (wa) obtained (a) from TB 
calculations with t — 2.70 (eV) and (b) from first-principles 
calculations (symbols). The solid lines in (b) are from Eq. 
(1). (c) First-principles band structures of A^a-AGNRs with 
Na~12, 13, and 14, respectively. 

exhibit direct bandgaps at kda = for all cases. 

A determining factor in the semiconducting behavior of 
A'q-AGNR is quantum confinement which can be charac- 
terized by Aat^ w-'^ 0,ll,[i,0,S,[i,[l3- In addition, 
as will be discussed below, the edge effects play a cru- 
cial role and force the (3p -I- 2)-AGNRs (predicted to be 
metallic by TB model) to be semiconductors. The edge 
carbon atoms of our AGNRs are passivated by hydro- 
gen atoms (by some foreign atoms/molecules in general) 
so that the a bonds between hydrogen and carbon and 
the onsitc energies of the carbons at the edges would be 
different from those in the middle of the AGNRs. The 
bonding distances between carbon atoms at the edges 
are also expected to change accordingly. Such effects 
have been observed in large aromatic molecules, e.g., ova- 
lene (C32H14) [2^. In Fig. 3(a), we show that the bond 
lengths parallel to dimer lines at edges (ai and a^^ for 
A^a-AGNR in Fig. 1 (a)) are shortened by 3.3-3.5% for 
the 12-, 13-, and 14-AGNR as compared to those in the 
middle of the ribbon. From the analytic expressions for 
TB matrix elements between carbon atoms in Ref. 27, a 
3.5% decrease in interatomic distance from 1.422 A would 
induce a 12% increase in the hopping integral between tf- 
orbitals. 

To see the consequence of such effects more clearly, we 
introduce a lattice model [Fig. 3(b)] which is equivalent 
to the AGNRs within the TB approximation 0,[i,[i,01- 
The set of eigenvalues of a brick type lattice shown in 
Fig 3.(b) at kda = is further equivalent to that of 
a two- leg ladder system with Na rungs [E IE 01- 
The Hamiltonian of this simpler model can be written 

as H = Y.n=l £^.,nal^na^,,n - Y.n=l tni<A,nO'2,n + 

^■c-) -J2n=i^ J2l=itn,n+iiai,naf.,n+i + i^-c.), where {n, 

fi} denote a site, e^,n site energies, t|| ^^-^ and t:^ the 
nearest neighbor hopping integrals within each leg and 
between the legs respectively, and a^_„ the annihilation 
operator of 7r-electrons on the n-th site of the /z-th leg. As 
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FIG. 3: (a) The ratio of the calculated change in the carbon- 
carbon distance ((5a„) (see Fig. 1(a)) to the carbon-carbon 
distance in the middle of the TVa-AGNRs, i.e., Sa^ = 100 x 
""""^ where Oc = ae = ay = 1.424 A for Na = 12, Uc = 
ar = 1.422 A for Na = 13, and = ar = as = 1.423 A for 
Na = 14, respectively, (b) Topologically equivalent structure 
to the Ai'a-AGNR shown in Fig. 1(a). For the especial case of 
fc = 0, a lattice with periodic ladders (left) can be folded into 
a two- leg ladder with A'^a rungs (right). 

discussed above and shown in Fig. 3(a), and e^.„ at 
the edges would differ from those in the middle of GNRs. 
Hence, considering the simplest but essential variation 
from the exact solvable model to approximate the realis- 
tic situations, we assume that tj^ — tj;^ = {1 + 6)t and 

t:^ = t for n — 2, ■ ■ ■ , Na — 1, and t\[ ^^^^ = t for all n's. 
The site energies are set at e^^„ = eg for n = 1 and Na 
and otherwise regardless of fj,. This model Hamiltonian 
is solved pertubatively and the resulting energy gaps to 
the first order in S and Sq are as follows, 
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where Ag^, A3p_|_;^ and A3p_|_2 are the gaps of the ideally 
terminated ribbon when S — Eq = 0. They are given by 
t 2 — 4 cos ^o'^il'^ and respectively. 
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The zeroth-order gaps are identical to the values obtained 
from numerical calculations in Fig. 2(a) 0,11, SQ- With 
t = 2.7 (eV) [Hi and i5 = 0.12, the calculated gaps ob- 
tained using Eq. (1) are in good agreement with our LDA 
results [Fig. 2(b)]. This implies that the 12% increase of 
the hopping integrals between carbon atoms at the edges 
opens the gaps of the {3p + 2)-AGNRs and decreases (in- 
creases) the gaps of 3p-AGNRs ((3p + l)-AGNRs). This 
analysis provide the physical explanation of the changes 
in the gap hierarchy discussed before. We note that there 
is no contribution from the variation in the site energies 
(eq) at the edges to first order. 

Next, we find that nanoribbons with zigzag shaped 
edges also have direct bandgaps which decrease with in- 
creasing width {wz). The eigenstates of the ZGNRs near 
Ep , without considering spins, have a peculiar edge-state 
structure. As noted earlier within the tight-binding pic- 
ture 0], there are two edge states decaying into the cen- 
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FIG. 4: (color online) (a) Contour graph for pa{r) ~ p^ir) of a 
12-ZGNR (the density is integrated over the normal direction 
to the ribbon plane). The lowest (highest) contour of ±0.4 x 
10~^aQ ^ is drawn by a thick blue (red) line and the spacing for 
blue (red) lines is 1.0 x 10~*aQ ^ (ao = Bohr radius), (b) The 
band structure of a 12-ZGNR. The ct- and /9-spin states are 
degenerate in all energy bands. and A^ denote the direct 
bandgap and the energy splitting at kdz = tt, respectively, 
(c) The variation of A° and A?, as function of the width (w^) 
of AT^-ZGNRs. The solid hue is a fit curve for A° {N^ > 8) 
and the dotted lines are drawn to guide the eyes. 

ter of the ZGNR with a decay profile depending on their 



momentum as e where ak 



s/3d. 
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< kdz < T^, dz — unit cell length shown in Fig. 1(b)). 
Our first-principles calculation also predicts a set of dou- 
bly degenerate flat edge-state bands at Ep when not con- 
sidering spins (not shown here). Since the edge-states 
around Ep form flat bands, they give rise to a very large 
density of states at Ep. Thus infinitesimally small onsite 
repulsions could make the ZGNRs magnetic unlike 
the case with two dimensional graphene which has a zero 
density of states at Ep. As pointed out in a TB study 
earlier 14 and later confirmed by first-principles stud- 
ies [III [3 ) our LSDA calculation also shows that the 
ground state of ZGNRs with hydrogen passivated zigzag 
edges indeed have finite magnetic moments on each edge 
with negligible change in atomic structure 11, 13, 13, l3 |. 

Upon inclusion of the spin degrees of freedom within 
LSDA, the ZGNR are predicted to have a magnetic insu- 
lating ground state with ferromagnetic ordering at each 
zigzag edge and antiparallel spin orientation between the 
two edges. The spatial spin distributions of the ground 
state in the case of 12-ZGNR is displayed in Fig. 4(a). 
The small spin-orbit coupling [i^ in carbon atoms is 
neglected in the present study, and we label one spin 
orientation as a-spin (red) and the opposite as /3-spin 
(blue) in Fig. 4 (a). The total energy difference per 
edge atom between nonspin-polarized and spin-polarized 
edge states increases from 20 meV (Nz = 8) to 24 meV 
{Nz = 16). These energy differences are further sta- 
bilized by an antiferromagnetic coupling between the 
two edges. The total energy difference between ferro- 
magnetic and antiferromagnetic couplings between edges, 
however, decreases as Nz increases and eventually be- 
comes negligible if the width is significantly larger than 
the decay length of the spin polarized edge states [l^ . 
The ferromagnetic-antiferromagnetic energy differences 
per unitcell are 4.0, 1.8, and 0.4 meV for the 8-, 16-, and 
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32-ZGNR, respectively. Our LSDA results agree with 
previous studies 0, 11, [l3 and consistent with a theo- 
rem based on the Hubbard Hamiltonian on a bipartite 
lattice poj . Though infinite range spontaneous magnetic 
ordering in a one dimensional Heisenberg model is ruled 
out [s^l , spin orderings in nanoscale system are realizable 
in practice 31, 32, 33, 34 1 at finite temperature assisted 
by the enhanced anisotropy on substrates 3^, 34 1 . 

The energy gaps in ZGNRs originate from the stag- 
gered sublattice potentials resulting from the magnetic 
ordering, which introduce bandgaps for electrons on a 
honeycomb hexagonal lattice [35|. This is realized be- 
cause the opposite spin states on opposite edges occupy 
different sublattices respectively (black rectangles on the 
left side and grey ones on the right belong to different 
sublattice respectively in Fig. 1(b)). So, the ZGNRs can 
be considered as the magnetic analog of a single BN sheet 
because the former has a bandgap which originates from 
the exchange potential difference on the two sublattices 
while the bandgap of the latter is from the ionic poten- 
tial difference between boron and nitrogen atoms located 
on the different sublattices The Hamiltonian (H) 

and Bloch wavefunctions {ipnka{i3)) of ths ground states 
satisfy [Ti,TM] — and TMipnka = ipnkp where T is 
the time-reversal symmetry operator and M. a mirror 
symmetry operator interchanging sites on opposite sides. 
Hence, a- and /3-spin states are degenerate in all bands 
and have the same gap as shown in Fig. 4(b). 

Since the strength of the staggered potentials in the 
middle of the ribbon decreases as the width increases, the 
LSDA bandgaps of ZGNRs are inversely proportional to 
the width. The calculated energy gaps (in eV) can be 
fit by Al{wz) = 9.33/(w^ + 15.0) with Wz in Angstroms 
when iV^ > 8 as shown in Fig. 4(c). It is also shown that 
the energy splitting at fca = tt (A?, in Fig. 4(b)) converges 
to 0.53 eV (32-ZGNR) from 0.51 eV (8-ZGNR). 

In summary, we have shown that graphene nanorib- 
bons with homogeneous armchair or zigzag shaped edges 
all have energy gaps which decrease as the widths of the 
system increase [37| . The role of the edges are crucial for 
determining the values and scaling rule for the bandgaps. 
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